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Our aim is to identify and classify mean-motion resonances (MMRs) for the coplanar circular restricted three-body prob- 
lem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which 
is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what 
kind of resonances are present, and the inspection of the orbital elements to classify the periodicity. We show that the 
2:1 resonance occurs the most frequently. Among other resonances, the 3:1 resonance is the second most common, and 
furthermore both 3:2 and 5:3 resonances occur more often than the 4:1 resonance. Moreover, the resonances in the copla- 
nar CR3BP are classified based on the behaviour of the orbits. We show that orbital stability is ensured for high values 
of resonance (i.e., high ratios) where only a single resonance is present. The resonances attained are consistent with the 
previously established resonances for the Solar System, i.e., specifically, in regards to the asteroid belt. Previous work 
employed digital filtering and Lyapunov characteristic exponents to determine stochasticity of the eccentricity, which is 
found to be consistent with our usage of Lyapunov exponents as an alternate approach based on varying the mass ratio 
instead of the eccentricity. Our results are expected to be of principal interest to future studies, including augmentations to 
observed or proposed resonances, of extra-solar planets in binary stellar systems. 
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1 Introduction 

The circular restricted three-body problem (CR3BP) has been 
considered in orbital mechanics for more than two centuries, 
and its numerous applications to describe the motion of plan- 
ets, asteroids and artificial satellites in the Solar System are 
well-known dDyorald|l98ill986HMurrav & Dermotfe oOO: 



Rabl & Dvorak|[l98a ISmith & Szebehelvll 19931 ISzebehelv 



1967 and references therein). More recently, the study of the 



CR3BP has been motivated as a testbed for investigations of 
the orbital stability of extra-solar planetary systems around 
single stars as well as planets in stellar binary systems. 

Extensive studies of orbital stability in the CR3BP have 
been performed by many authors , who used tools ranging 
from the Jacobi constant (e.g.. 



Cunt z. Eberle & Musiela 



2008), a hodograph-based 



The origin and nature of mean-motion resonances (MMRs) 
in the CR3BP, particularly for cases of low mass ratios, have 
also been studied, and the role of MMRs in orbital stabil- 
ity has been investigated in detail. iHadjidemetriou dl993 ) 
identified three (2:1, 3:1 and 4:1) main MMRs in the Sun- 
Jupiter and asteroid system. Extensive studies of MMRs, 

with extensions to the elliptic rest ricted three-body pro b lem 

(ER3BP), have been per f ormed bvlFerraz-Mellol d 1 9941) ; lHaghi ghipour et al.1 
i2003[) :| Mardling|(l2007b:lNesvorny et al.ld2002h;|pilat-Lohinger & Dvorak! 
(200a: ISzenkovits & Makdld2008l) ; IVela -Arevalo & Marsden 
(|2004J). It has been shown that the main difference in orbital 
stability between stable and unstable systems with MMRs is 
that the latter involves more than one resonance, and that the 

existence of two or more resonanc es leads to their overlap 

makin g some orbits unstable (e.g.. lMardling|2007tlMudrvk & Wul 
2006). 



2007; Eberle, Cuntz & Musielak 

method (|Eberle & CjmfabQlOj) and Lyapunov exponents (e.g., of work is a i so motivated by earlier r esults on 

Gonczi & Froeschle 1981; Jefferys & ¥^1983; , Lecar. Frankhn.&Munso ^.^^ sta r _ plan et systems . Previously. tM^cTitlT 
1992;' : " ^ : — o 



2001 



Lissauerll999 



Quarl es et al 



, Milani & No bUai992tlMurrav & HolmarJ - 2001) and fi~ 
20111 [Smith & Szebehely|l993j). There ' 



have also been inquiries into different numerical methods 
speci fically designed to study the orbital stability problem 



et al.l (|2006) discovered two planets in 2: 1 
orbital resonance about the star GJ 876 and HD 82943, re- 
spectively. Other interesting cases include the 3:1 MMR in 



spci mccmy ucMyicu tu Muuy uic uiuiuu ^aumiy uiuuicm ^ fiye p j anet syfJtem Qf 55 Cnc Qccmling between the in- 

(e.gJDavid et al. 2003; Holman & WiegerrJ1999HMusielak et al^ planets 55 Cnc c and 55 Cn c b djj e t al.l2003HMarcv et al. 



120051; | Yoshida| | 1990| >. Some recent studies involve space- EooaiNovak, Lai & Linl2003l) and the 3:2 MMR in the sys 



craft trajectories to the Lagrangian point L 3 dTantardini et al. | tem f HD 45364 dCorreia et all EoOOh . Furthermore, the 



2010) and the stabilit y analysis of artificial e quilibrium points 
in the CR3BP l lBombardelli & Pelaez|[2oTl . 



Corresponding author: e-mail: cuntz@uta.edu 



two ou ter planets of ' 47 U Ma are found to be in a 5:2 reso- 
nance dFischer et al. 2002 ). This system is of particular in- 
terest as it is hosting two gas planets with a similar mass ra- 
tio than that between Jupiter and Saturn, which are orbiting 
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a star of similar spectra type than the Sun. However, the gas 
planets in this system are relatively close to the zone of hab- 
itability; therefore, they ad versely affect the astrobiological 



significance of the system (Cuntz et al. 2003c iGozdziewski 
20021: iNoble. Musielak & Cuntzj|2002h . 



To identify MMRs in the CR3BR particularly for small 
mass ratios, or in the elliptical restricted three-body prob- 
lem ( ER3BP) several differe nt methods have been devel- 
oped. Hadiide metrioul ( 1993 ) used the method of averaging 
as well as the method of computing periodic orbits numeri- 
cally. He showed how to utilize these methods to determine 
MMRs; furthermore, he discussed the advantages and dis- 
advantages of each method as well as the effects of their in- 
herent limitations on the attained results. Another method of 
finding MMRs is base d on the concept of Lyapunov expo- 
nent as introduced by INesvornv & Morbidellil dl998 ) with 



2 Theorectical Approach 

We consider the CR3BP with objects of masses mi, 777,2 and 
7713; note that mass 7773 is assumed to be negligible com- 
pared to mi and m 2 . From a mathematical point of view, 
the problem is de scribed by the following set of equations 
dSzebehelvlll967h 

.. „. 9fl .. „. on .. on 

x-2y=—, y + 2x = —, z = — , (1) 
ox ay oz 

where the potential function il is given by 



1999) 



its further development by iMorbidelli & Nesvornv 
and INesvornv et al.l (120021) . So far, the developed methods 
have mainly been applied to identify MMRs that affect the 
structure and evolution of the asteroid belt. This means that 
the identified MMRs were found for the specific mass ratios 
of the Solar System primaries (Sun and Jupiter, and Jupiter 
and Saturn). However, there is significant ongoing work as 
well as the need for future work pertaining to extra-solar 
planetary systems, which is self-evident from the recent dis- 
coveries by the Kepler space telescope of the circumb inary 
planets Kepler- 1 6 , Kepler-34 and Kepler-35 dPoyle et alJ201 1 
Welsh et alJEoil . 



The main purpose of this paper is to identify and clas- 
sify MMRs in the CR3BP for a broad range of mass ratios of 
the stellar components. Our method of finding resonances is 
based on the maximum Lyapunov exponent, which is used 
as an indicator for the location of resonances once the masses 
of the two massive bodies are specified. Then, a Fast Fourier 
Transform (FFT) is used to determine what kind of reso- 
nances exist. Thereafter, periodicity is determined through 
the inspection of the orbit diagrams and osculating semi- 
major axis of representative cases. Future applications of 
enhanced versions of our results will involve (1) extra-solar 
planetary systems around single stars (2) exomoons in star- 
planet systems, and (3) planets in stellar binary systems. 
Currently, the ideal application of the CR3BP to planets in 
stellar binary systems remains undetected. Therefore, spe- 
cific applications are beyond the scope of this paper that 
is fully devoted to the problem of finding and classifying 
MMRs in the CR3BP. Nonetheless, our present study is an 
important step toward envisioned future investigations mo- 
tivated by existing and anticipated discoveries of planets 
through ongoing and future search missions. 

The paper is structured as follows: Our method is de- 
scribed in Sect. 2 and in Sect. 3 we present our discus- 
sion and results. Our summary and conclusions are given 
in Sect. 4. 



ft = — + — + £ Ki + pr 2 2 ) 
n r 2 2 



1 



with fi = 777^/(7774 + m 'i)t ol 

n2 



1 — /j, and 

[x — p) 2 + y 2 + z 2 

2 1 1 \2 . 2 1 2 

r 2 = (x + a) +y + z . 
The Jacobian integral of Eq. (Q~|i is 

C = 2fl 



(x 2 +y 2 + z 2 ) 



(2) 



(3) 



(4) 



where C is the so-called Jacobi constant; see also previ- 
ous papers of this series for additional explanation on the 
adopted method. 

We study the manifestations of resonances in numeri- 
cal solutions in the physical plane (x, y) for the Copenhagen 
case (p = 0.5) and the special cases p = 0.1, 0.2, 0.3, and 
0.4. We vary the parameter po, which is defined by the ra- 
tio of the initial distance of mass 7713 from the host mass 
mi, given as Ro, and the separation distance, D, between 
the two larger (mi and m 2 ) masses, i.e., po = Rq/D. As 
an initial condition, the object of mass 777,3 starts at the 3 
o'clock position with respect to the massive host star, mi. 
We also assume a circular initial velocity for 7773. The nu- 
merical simulations are performed using a sixth order sym- 
plectic integration metho d with a consta nt time step of 10~ 3 
binary periods per step (lYoshidal 1 1 9901) . Unstable systems 
are confirmed using a sixteenth order Gragg-Burlisch-Stoer 
integration method with the aforementioned constant time 
step as an initial guess JGrazier et al. 19961). Add itionally, 
we use fitting formulas (Hol man & Wiegertlll999l) to con- 
firm the numerically determined stability limits. 

In order to determine resonances, we utilize the method 



of Lyapunov exponen ts (e.g., Nesvornv & Morbidellill998l 
Nesvornv et al.l 12002b . In this method, we consider a con- 



stant mass ratio p, calculate the spectrum of Lyapunov ex- 
ponents, and monitor the maximum Lyapunov exponent with 
respect to the parameter po . Although we calculate the full 
spectrum of Lyapunov exponents, we limit this discussion 
to the maximum Lyapu nov exponent for the reasons ex- 
plained elsewhere (e.g JOuarles et al l201lh . The maximum 
Lyapunov exponent is used as an indicator of the maxi- 
mum role that a particular resonance will have on the sys- 
tem. This provides us with a way of inspecting how the 
maximum Lyapunov exponent changes in response to small 
variations in the initial condition po- Addit ionally, we use 
the fast Lyapunov indicator (e.g jFroeschle. Gonczi & Legal 
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Fig. 1: Maximum Lyapunov exponent is plotted versus the distance ratio po for jj, = 0.1. The stable and unstable cases are 
denoted. 
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Fig. 2: Same as Fig. 1 but for // = 0.2. 



1997l:lLega & Froeschlel2 001) to estimate the extent of weak 
chaotic or ordered motion. 

After choosing a value of /i and performing calculations 
of the maximum Lyapunov exponent, we use a FFT to pro- 
duce a periodogram, and then take ratios of the periods in 
the periodogram to determine what kind of resonances (if 
any) are present. Using the Lyapunov exponent with a peri- 
odogram provides a method of quickly identifying regions 
of possible resonances rather than invoking a brute force 
method of performing a FFT over all possible values of po- 
The ratios of the periods in the periodogram are taken with 



respect to the period axis; however, the strength of the peaks 
are not considered. In the tables, we provide the periods of 
the three strongest peaks in the periodogram and label them 
as Peak 1, Peak 2, and Peak 3 in ascending order. Reso- 
nance 1 denotes the closest rational expression of the ra- 
tio of Peak 2 to Peak 1, whereas Resonance 2 denotes the 
closest rational expression of Peak 3 to Peak 1 . More peaks 
may be possible in the periodograms; however, the induced 
instability greatly overwhelms the necessity for investigat- 
ing them. There are also elements in the tables without data 
(shown by ellipsis). They indicate simulations where a Peak 
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Distance Ratio p Q 

Fig. 3: Same as Fig. 1 but for p = 0.3. 



3 could not be discriminated against the background; there- 
fore, a value for Resonance 2 could not be determined. We 
study the range of po beginning at po = 0.250 and ending 
at po = 0.700 in increments of 0.001 for a given value of p. 

3 Results and Discussion 

3.1 Case studies 

We consider five different cases of fixed mass ratio p, start- 
ing at p = 0.1 and proceeding in increments of 0.1. Each 
case presents a different aspect of the problem but together 
the combined results allow us to identify and classify res- 
onances in the coplanar CR3BP for relatively large values 
of p. For each case, we start with the parameter po at the 
value 0.250 that is increased in increments of 0.001 up to a 
final value of po = 0.700. We show one plot and one table 
for each case of p. In each plot, we distinguish between the 
stable cases, which continued for the full 10 5 binary orbits 
of simulation time by a solid line, and the unstable cases by 
dashed lines; the latter ended early due to the ejection of the 
planet from the system. Our criteria for ending each simu- 
lation are based on checking the error in the Jacobi constant 
and the value of the maximum Lyapunov exponent where 
high values of either im ply a high prob ability for to be 
ejected from the system (Ou arles et al.ll201 II) . We also pur- 
sued numerical comparisons with an independent integra- 
tion method to verify the integrity of the proposed outcome 
(see Section [2]). These separate criteria provide us with a 
method for checking for instability and stability, as well as 
a measure for identifying chaos in the system. In the ac- 
companying tables, we provide the values of po where the 
maximum Lyapunov exponent is found to peak. 

After the maximum Lyapunov exponent is calculated, 
we transform the system pointwise from barycentric coor- 



dinates to Jacobi coordinates choosing the more massive 
body as reference point. We perform a FFT using the time 
series data for that specific po- With this result, we obtain 
the Fourier spectrum and convert it to a corresponding peri- 
odogram. The resonances are determined by taking the ratio 
of periods for the strongest peaks. In the tables, we provide 
the computational results for the three strongest peaks, if 
available, for that value of po. The resonances are given as 
a ratio of the period of the small (m^) mass to the period of 
the smaller binary (m2) mass. For example, 4:1 would tell 
us that the small mass orbits four times in a single orbit of 
the smaller binary mass. 

The first case of p = 0.1 is depicted in Fig. Q] with nu- 
merical results given in Table Q] This case illustrates a base- 
line of configurations where the second large (TO2) mass is 
not greatly perturbing the small (m^) mass from po = 0.250 
to 0.350. This provides a clear indication of where a possi- 
ble resonance is not substantial enough to alter the orbit. 
However, the value of the maximum Lyapunov exponent 
peaks first at po = 0.394. This indicates that the second 
large mass is at a specific geometric position allowing it to 
greatly perturb the motion of the small mass. The same fea- 
ture occurs again at po = 0.457 and 0.471. Note that the 
resonances associated with each peak in Fig. Q]are n °t all 
the same; see Table [T] Moreover, we identify a large region 
of instability when po is further in creased, which is con- 
sistent with previou s investigations ( Eberle & Cuntz 2010; 
Ouarles etai]l201lb . 

The second case of p = 0.2 is depicted in Fig. |2] and 
the relevant numerical values are presented in Table This 
case differs from the preceding one as now the first reso- 
nance peak occurs sooner at po = 0.334. This is expected 
because by increasing the value of p, the perturbing mass 
has a greater ability to interact with the smaller mass to al- 
ter its orbit. There exists a stability region from po = 0.250 
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Table 1: Resonances for selected values of /?o with an elevated maximum Lyapunov exponent given for p = 0.1. The 
different orbital types are classified as Periodic, Quasi-Periodic, and Non-Periodic, and are denoted by P, QP, and NP, 
respectively. 



Po 


Peak 1 


Peak 2 


Peak3 a 


Resonance 1 


Resonance 2° 


Orbital Type 


0.394 


0.3494 


1.045 




3:1 




P 


0.457 


0.4693 


0.8827 


1.001 


15:8 


17:8 


P 


0.471 


0.4832 


0.9671 




2:1 




P 



"Note: Elements without data indicate simulations that did not have a Peak 3 in the periodogram; therefore, a value for Resonance 
2 could not be determined. 







Table 2: Same as Table 1 but for 


p — U.Z. 




PO 


Peak 1 




Peak 3 


Resonance 1 


Resonance 2 


Orbital Type 


0.334 


0.2715 


1.087 




4:1 




P 


0.391 


0.3631 


1.143 




22:7 




QP 


0.399 


0.3756 


1.127 




3:1 




QP 


0.450 


0.4636 


0.7716 


1.162 


5:3 


5:2 


QP 


0.458 


0.4769 


0.815 


1.151 


12:7 


12:5 


OP 


0.482 


0.5153 


1.031 




2:1 




P 


0.491 


0.3457 


0.5186 


1.037 


3:2 


3:1 


QP 






Table 3: 


Same as Table 1 but for 


p — U.O. 




PO 


Peak 1 


Peak 2 


Peak 3 


Resonance 1 


Resonance 2 


Orbital Type 


0.355 


0.3267 


1.297 




4:1 




QP 


0.370 


0.3503 


1.357 




31:8 




QP 


0.454 


0.4743 


1.448 




3:1 




QP 


0.471 


0.4996 


1.500 




3:1 




P 


0.483 


0.5121 


0.7680 


1.537 


3:2 


3:1 


OP 


0.488 


0.5244 


0.8034 


1.511 


3:2 


23:8 


NP 


0.493 


0.5236 


0.8278 


1.495 


14:3 


14:5 


QP 




u.jyuo 


0.9875 


1.472 


5:3 


J.Z 


QF 


0.537 


0.6016 


1.017 


1.475 


5:3 


22:9 


QP 






Table 4: Same as Table 1 but for 


p = 0.4. 




Po 


Peak 1 


Peak 2 


Peak 3 


Resonance 1 


Resonance 2 


Orbital Type 


0.429 


0.4456 


1.995 




9:2 




QP 


0.444 


0.4726 


1.732 




11:3 




QP 


0.471 


0.5236 


1.580 




3:1 




QP 


0.505 


0.2958 


0.5918 


1.480 


2:1 


5:1 


QP 


0.507 


0.2976 


0.5953 


1.473 


2:1 


5:1 


QP 


0.518 


0.2046 


0.3070 


0.614 


3:2 


3:1 


QP 


0.521 


0.2060 


0.3090 


0.618 


3:2 


3:1 


QP 






Table 5: 


Same as Table 1 but for 


p = 0.5. 




Po 


Peak 1 


Peak 2 


Peak 3 


Resonance 1 


Resonance 2 


Orbital Type 


0.410 


0.232 


0.4639 




2:1 




P 


0.426 


0.251 


0.5014 


1.506 


2:1 


6:1 


QP 
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Fig. 4: Same as Fig. 1 but for p = 0.4. 
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Fig. 5: Same as Fig. 1 but for fi = 0.5. 



to 0.333. However, the value of the maximum Lyapunov ex- 
ponent peaks at seven different values of pg, which demon- 
strates the greater influence of the perturbing mass m 2 com- 
pared to the case of p = 0.1. Also, there are apparent gaps 
surrounding the instability regions. Actually, they are due to 
our increment size in po and, therefore, these regions should 
be treated as continuous with appropriate lines transition- 
ing from one region to the next. There is also a resonance 
plateau in the range of po — 0.45 to 0.48. This illustrates 
a region where the phenomena of resonance is widespread; 
note there is a similar feature in the study of asteroids (i.e., 
the 3:1 resonance at 2.5 AU or po = 0.48). 

The third case of p — 0.3 is depicted in Fig. [3] and the 
associated numerical results are given in Table [3] Here the 



first resonance peak occurs at po = 0.355. It also exhibits a 
well-pronounced region where the resonance is not substan- 
tial enough to alter the orbit of the small mass 7713 as shown 
in Fig.[3]from po = 0.250 to 0.354. The value of the max- 
imum Lyapunov exponent peaks at 9 different values of po 
be fore reaching the stability limit at po = 0.600 as shown 
bv lEberle et alj d2008l) . It is important to note that the num- 
ber of resonance peaks has increased thus far as the value 
of p has increased. Also this case includes a peak in the 
maximum Lyapunov exponent, where the system is unstable 
at po = 0.488 leading to planetary ejection. The transition 
to this instability is shown by the developing resonances at 
po = 0.483 in Tabled 
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The fourth case of p = 0.4 is depicted in Fig. [4] with 
numerical results given in Table |4] This case differs from 
the preceding trend by the first resonance peak; this peak is 
delayed and occurs at po = 0.429 after an instability island. 
This is due to the perturbing mass gaining an overwhelming 
amount of gravitational influence over the small mass and 
creating an instability region from po — 0.350 to 0.400. 
Due to this change there are fewer positions for resonance 
to occur, the resonances become denser with respect to p , 
and the value of the maximum Lyapunov exponent peaks at 
7 different values of po. 

The final case of p = 0.5 is depicted in Fig. [5] and the 
relevant numerical results are given in Table |5] This case 
differs as the first resonance peak is found to occur sooner 
at po = 0.410. The Copenhagen case shows the worst case 
scenario for the small mass because the tug of war between 
the large masses becomes the most extreme. As a result 
there is only two narrow regions of stability. The value of 
the maximum Lyapunov exponent peaks at only two differ- 
ent values of pq. The first resonance characterizes a stable 
resonance, whereas the second is a resonance that shows the 
transition from stability to instability via resonance overlap. 

3.2 Classification of resonances 

In Tables 1 to 5, we identified the existence of 2 possible 
resonances. We now discuss them separately starting with 
the results obtained for Resonance 1. These results show 
that the 2:1 resonance occurs in all considered cases of p, 
except p = 0.3. It is also shown that the 2: 1 resonance is the 
most common resonance in the coplanar CR3BP in agree- 
ment with previous asteroidal investigations. The 3:1 reso- 
nance is the second most common resonance as it occurs in 
all cases except p = 0.5. A new and interesting result is that 
both the 3:2 and 5:3 resonances occur more often than the 
4:1 resonance, which only occurs for p = 0.2 and 0.3. 

The overall picture obtained from the special cases of 
the values considered for p is that the best chance for a res- 
onance occurs at the intermediate value of p = 0.3. In this 
case, we see that systems start with a resonance reflecting a 
large ratio of the small mass period to the large mass period. 
Then as the parameter po increases, so does the value of the 
ratio. After crossing the smoothest peak at po = 0.454, the 
secondary resonances are formed; these resonances are la- 
beled as Resonance 2 in Tables 1 to 5 and once they appear 
their existence may imply instability caused by resonance 
overlap. Among these resonances, the 3:1 resonance is the 
most dominant but new resonances such as 5:2 also occur, 
which is present for p = 0.2 and 0.3, and 5:1 occurring 
for p — 0.4. An interesting result is that only the 3:1 res- 
onance occurs also as a secondary resonance and that none 
of the other resonances labeled as Resonance 1 also occur 
as Resonance 2. 

Additionally, we can classify resonances based on the 
behaviour of the orbits. In Tables Q] to [5] we show classifi- 
cations based on orbit diagrams in rotating coordinates (see 



Figs.|6]to|9^) and the osculating semimajor axis of the small 
mass (see Figs. [6] to [9})). The three classes are periodic, 
quasi-periodic, and non-periodic, denoted by P, QP, and NP, 
respectively. Periodic orbits are those that are well-defined 
(i.e., minimal deviation from the path traced by previous 
orbits). They are closed, show periodic variations in the os- 
culating semimajor axis, and show two dominant peaks in 
the corresponding periodogram (see Figs. [5J;, 0;, and |9j;). 
Quasi-periodic orbits are those that precess or fill an an- 
nulus with an inner and outer boundary. Moreover, quasi- 
periodic orbits exhibit a more complicated osculating semi- 
major axis as well as a periodogram with several peaks (see 
Figs. |3),c). Non-periodic orbits are those that do not show 
a regular pattern and are thus most likely to be unstable. 
Also non-periodic cases show the greatest number of peaks 
in their periodograms. Of the cases presented, a low value 
of p yields a greater probability for periodic orbits. As the 
mass ratio p increases, the probability for quasi-periodic or- 
bits increases due to increasing influence of the perturbative 
mass and duration of the perturbation. However, there are 
some periodic orbits that exist for large values of p. They 
occur near the smooth plateau that forms in each of the fig- 
ures. 



4 Summary and conclusions 

The case studies of the coplanar CR3BP presented in this 
paper were used to identify and classify resonances in the 
coplanar CR3BP. Our method of finding resonances is based 
on the concept of Lyapunov exponents, or more specifically, 
it uses the maximum Lyapunov exponent. Based on the lat- 
ter, we were able to determine where the resonances oc- 
curred for a fixed value of mass ratio p. The conditions for 
uncovering these resonances were determined by the value 
of the maximum Lyapunov exponent at the end of each sim- 
ulation. The fact that the value of the maximum Lyapunov 
exponent peaks at specific values of distance ratio po pro- 
vides a strong indicator for the positions of the resonances. 
Using this method, we have shown that stability is ensured 
for high values of resonance (i.e., high integer ratios) where 
only a single resonance is present. 

We identified the primary and secondary resonances in 
the coplanar CR3BP as Resonance 1 and 2 (see Tables 1 
to 5), have shown that the 2:1 resonance occurs most fre- 
quently as Resonance 1 in the considered systems, and have 
shown that the 3:1 resonance is the next most common. 
An interesting result is that there are two other resonances, 
which are 3:2 and 5:3, and that they occur more often than 
the 4:1 resonance. Among the resonances labeled as Reso- 
nance 2, the 3:1 resonance is the most dominant but there 
are also new resonances such as 5:2 and 5:1. Moreover, we 
found that the 3:1 resonance is the only one that also occurs 
as a secondary resonance, and that none of the other reso- 
nances identified by us as Resonance 1 also occur as Res- 
onance 2. The existence of the resonances labeled as Reso- 
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Fig. 6: Case study showing the results for p = 0.2 and p = 0.334. Panel (a) shows the orbit of a planet in a rotating 
coordinates system, (b) shows the osculating semimajor axis for the first 50 binary orbits, and (c) shows the Fourier 
periodogram to determine the possible resonances. 
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Fig. 7: Case study showing the results for p = 0.3 and p = 0.471. Panel (a) shows the orbit of a planet in a rotating 
coordinates system, (b) shows the osculating semimajor axis for the first 50 binary orbits, and (c) shows the Fourier 
periodogram to determine the possible resonances. 



nance 2 is important as it may imply the onset of instability 
caused by the resonance overlap. 

We also classified the resonances in the coplanar CR3BP 
based on the behaviour of the orbits. The following three 
different classes were used : periodic (P), quasi-peri odic (QP), 
and no n-periodic (NP); see Eberle & Cuntz ( 2010l) and Quarles 
(1201 II) for previous usages of these classes in conjunction 
with proposed criteria for the onset of orbital instability. 
Among the considered cases, a low value for p and po en- 
tailed a greater probability for periodic orbits to occur. Then, 
the probability for quasi-periodic orbits increases as p in- 
creases. However, we also discovered some periodic orbits 
that exist for large values of p. 

The resonances reported in this paper are consistent with 
the previously established resonances for the Solar System, 
specifi cally, in the astero id belt with its dominant 2:1 reso- 
nance (Ferraz-MeloJll994|). They used digital filtering and 
Lyapunov characteristic exponents to determine stochastic- 



ity of the eccentricity. This is consistent with our usage of 
Lyapunov exponents as we took an alternate approach, which 
is varying the mass ratio instead of the eccentricity. Our ap- 
proach provides general results for the cases where p is con- 
siderably greater than the value of p in Sun-Jupiter-type 



s^s^Ems but where the eccentricity is small. Although our 
results have been obtained for the special case of the copla- 
nar CR3BP, we expect that it will be possible to augment our 
findings to planets in general stellar binary systems. Desired 
generalizations should include studies of the ER3BP (e.g. , 
Pilat-Lohinger & Dvorakll2002l ISzenkovits & Mak6ll2008l) 
as well as applications to planets in previously discovered 
star-planet systems. 
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(a) (b) (c) 

Fig. 8: Case study showing the results for p = 0.4 and p = 0.471. Panel (a) shows the orbit of a planet in a rotating 
coordinates system, (b) shows the osculating semimajor axis for the first 50 binary orbits, and (c) shows the Fourier 
periodogram to determine the possible resonances. 
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Fig. 9: Case study showing the results for p = 0.5 and p = 0.410. Panel (a) shows the orbit of a planet in a rotating 
coordinates system, (b) shows the osculating semimajor axis for the first 50 binary orbits, and (c) shows the Fourier 
periodogram to determine the possible resonances. 
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